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We use (3+l)-dimensional hydrodynamics with exact lon- 
gitudinal boost-invariance to study the influence of collision 
centrality and initial energy density on the transverse flow 
pattern and the angular distributions of particles emitted near 
midrapidity in ultrarelativistic heavy-ion collisions. We con- 
centrate on radial flow and the elliptic flow coefficient V2 as 
functions of the impact parameter and of the collision en- 
ergy. We demonstrate that the finally observed elliptic flow 
is established earlier in the collision than the observed radial 
flow and thus probes the equation of state at higher energy 
densities. We point out that a phase transition from hadronic 
matter to a color-deconfined quark-gluon plasma leads to non- 
monotonic behaviour in both beam energy and impact param- 
eter dependences which, if observed, can be used to identify 
such a phase transition. Our calculations span collision en- 
ergies from the Brookhaven AGS (Alternating Gradient Syn- 
chrotron) to beyond the LHC (Large Hadron Collider); the 
QGP phase transition signature is predicted between the low- 
est available SPS (CERN Super Proton Synchrotron) and the 
highest RHIC (Brookhaven Relativistic Heavy Ion Collider) 
energies. To optimize the chances for applicability of hydro- 
dynamics we suggest to study the excitation function of flow 
anisotropies in central uranium-uranium collisions in the side- 
on-side collision geometry. 



I. INTRODUCTION 

At a given beam energy, the highest energy densities 
can be reached in central collisions (impact parameter 
b = 0) between the largest available nuclei. Hence for 
many years the experimental and theoretical attention 
has focussed on such collisions. Non-central (b ^ 0) col- 
lisions are, however, interesting in their own right since 
they exhibit new phenomena which are forbidden by azi- 
muthal symmetry in central collisions between spherical 
nuclei. For non-central collisions the directions of the 
beam axis and the impact parameter b define the collision 
plane, and many interesting physical phenomena are now 
non-trivial functions of the azimuthal angle ip relative to 
the collision plane. These include in particular the trans- 
verse geometry of the collision fireball as measured with 
two-particle Bose-Einstein correlations (see e.g. Q and 
references therein) and momentum-space anisotropies in 
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the transverse plane due to anisotropic transverse flow of 
the fireball matter Q . 

Aside from changing the collision energy, limited vari- 
ations of the energy density of the reaction zone are also 
possible by varying the collision centrality. Variation of 
the initial energy density provides the handle for study- 
ing phase transitions in nuclear matter, in particular 
the quark-hadron transition at a critical energy density 
e c <lGeV/fm 3 ||. Non-central collisions between spher- 
ical nuclei and/or central collisions between deformed nu- 
clei provide new opportunities to correlate phenomena 
related to azimuthal anisotropies with the initial energy 
density. This may yield novel phase transition signatures. 
In Q this idea was exploited for the so-called directed 
flow at forward and backward rapidities: the softening of 
the equation of state (EOS) in the phase transition region 
was predicted to lead to a reduction of the directed flow, 
making the phase transition visible as a minimum in its 
excitation function. Sorge Mj^ suggested analogous fea- 
tures for the elliptic flow @-^| which were further studied 
in 1 10 12| . The effects of a phase transition on the exci- 
tation function of radial flow in central collisions between 
spherical nuclei had been discussed earlier in |l3|-|l5| . 

An important difference between the radial flow ob- 
served in azimuthally symmetric central collisions and 
the anisotropic directed and elliptic flows in non-central 
collisions and/or central collisions between deformed nu- 
clei was pointed out by Sorge in ||: 

1. Directed flow affects mostly particles at forward and 
backward rapidities which (at energies above a few hun- 
dred MeV/nucleon) are deflected away from the beam 
direction by the pressure built up between the colliding 
nuclei during the time of their mutual overlap. Since 
the thus affected particles quickly leave the central re- 
gion where this transverse pressure force acts, the finally 
observed directed transverse flow pattern is established 
very early in the collision. Its natural time scale is given 
by the transition time of the two colliding nuclei which 
decreases with increasing beam energy; this causes a de- 
crease at high collision energies (after an initial rise at 
low beam energies) of the directed flow ||] . This decrease 
is amplified by a lack of thermalization during the very 
earliest stages of the collision which prohibits fast enough 
buildup of transverse pressure and thus eventually invali- 
dates the applicability of hydrodynamic concepts for cal- 
culating the directed flow. Such pre-equilibrium features 
may even cover up Q the phase transition signal Q in 
the excitation function of directed flow. 

2. The elliptic flow is strongest near midrapidity [ jT7| . 
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Its driving force is the azimuthal anisotropy of the trans- 
verse pressure gradient, caused by the geometric defor- 
mation of the reaction region in the transverse plane. As 
pointed out in ||jll]], elliptic flow acts against its own 
cause by eliminating the geometric deformation which 
generates it, thereby shutting itself off after some time. 
This time is, at least at high energies, longer than the 
nuclear transition time. Elliptic flow is thus generated 
later than directed flow, and hydrodynamic concepts for 
its description may have a larger chance of being valid, 
even if the spatial deformation which causes elliptic flow 
exists only for a fraction of the total fireball lifetime. An 
important focus of this work will be a quantitative deter- 
mination of the time scale over which elliptic flow is gen- 
erated, as a function of the collision energy. We will see 
that this time scale grows with the overall size of the (ini- 
tially deformed) collision region Studying central 
collisions between large deformed nuclei like 238 U [ ^8|Jl9| 
therefore improves the chances that thermalization hap- 
pens sufficiently early for a hydrodynamic description of 
elliptic flow evolution to be valid. Such collisions are the 
preferred proving ground for hydrodynamic predictions 
for the excitation function of elliptic flow. 

3. Radial flow is generated by the pressure gradient 
between the interior of the collision fireball and the ex- 
ternal vacuum; this force persists throughout the fireball 
expansion until freeze-out. Of all three transverse flow 
patterns it thus has the strongest weight at late times. 
Comparing the excitation functions of elliptic and radial 
flow with their intrinsically different weights for the EOS 
at early and late times (i.e. at high and low energy den- 
sity) may help with the identification of phase transition 
signatures and their discrimination against possible non- 
equilibrium effects from incomplete local thermalization. 
Of course, the final proof for the phase transition to quark 
matter will require an additional correlation of the here 
predicted structures in the anisotropic flow pattern with 
other "quark-gluon plasma signatures" (see ) . 

As already indicated we here study the evolution of 
transverse flow in a macroscopic hydrodynamic frame- 
work (to be contrasted with microscopic kinetic ap- 
proaches [|]]|,[l(|). This approach, which is based on the 
assumption of rapid local thermalization, allows the most 
direct connection of observables to the EOS of the hot 
matter in the collision fireball, including possible phase 
transitions. Its validity can be tested both experimen- 
tally and by comparison with kinetic approaches. We will 
not do so here (see, for example, Rcfs. [ [10|]24| ) but rather 
concentrate on qualitative predictions resulting from the 
hydrodynamic approach. 

Hydrodynamics cannot describe the earliest collision 
stage of nuclear energy loss and entropy production by 
thermalization of the energy deposited in the reaction 
zone by the stopping process; this must be replaced by 
appropriate initial conditions for the hydrodynamic ex- 
pansion. The evolution of azimuthally asymmetric reac- 
tion zones requires a (3+l)-dimensional hydrodynamic 
approach. This is very time consuming and makes a tun- 



ing of initial conditions to data difficult [£5|,g6|]. However, 
near midrapidity and especially for high collision energies 
the longitudinal expansion dynamics is expected to be 
given by the Bjorken scaling solution |22| which can be 
implemented analytically. The remaining hydrodynamic 
equations for the transverse dynamics live in 2 space and 
1 time dimension and are much easier to solve [j|,|ll],|l2| . 
The hydrodynamic evolution is terminated by a freeze- 
out criterium (in our case a fixed decoupling energy den- 
sity). At this point the energy and baryon densities are 
converted into temperature and chemical potentials for 
baryon number and strangeness, using the EOS, and the 
particle spectra are calculated using the Cooper- Frye pre- 
scription j2^] . With these spectra and the hydrodynamic 
flow pattern on the freeze-out surface the average radial 
flow velocity ({v±)) and the elliptic flow coefficient V2 are 
evaluated. 

The present paper gives technical details for our pre- 
vious two short reports in flH| and significantly extends 
the results presented there. The excitation function for 
v 2 is complemented by a similar one for the average ra- 
dial flow and calculated up to very much higher energies. 
We also compute the impact parameter dependence at 
fixed beam energy of the elliptic flow scaled by the ini- 
tial spatial anisotropy. The time evolutions of radial and 
elliptic flow and their dependence on the collision energy 
are discussed in detail, in order to establish to which ex- 
tent elliptic flow is really a signature for early pressure 
in the system ]5|,^|,[ll| . Finally, we advertize central U+U 
collisions in the side-on-side configuration as an optimum 
system for studying the hydrodynamic evolution of ellip- 
tic flow and the quark-hadron phase transition signature 
in its beam energy dependence |ll ] . We give predictions 
for the time evolution of radial and elliptic flow, for their 
excitation function and for the p T -dependence of the el- 
liptic flow coefficient at SPS energies for this particular 
collision system. 



II. THE HYDRODYNAMIC MODEL 

The equations of relativistic ideal hydrodynamics fol- 
low from the local conservation laws for energy, momen- 
tum, and other conserved currents (e.g. baryon number), 

d^ u (x) = and V(s)=0, (2.1) 

by inserting the ideal fluid decompositions 

T^{x) = (e(x) + p{x]}u»{x)u v {x) - g^p(x) , (2.2) 



j^(x) = n(x) ^(x) 



(2.3) 



e(x) is the energy density, p{x) the pressure, and n{x) 
the conserved number density at point x^—{t, x, y, z); 

u' J, (x)=j(l,v x ,v y ,v z ) with 7= 1/ ^jl—v%—vl—vl is the 

local four velocity of the fluid. Ideal hydrodynamics as- 
sumes that local thermalization by the strong interac- 
tions among the matter constituents happens fast on the 
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scale defined by the space-time gradients of these quan- 
tities and therefore neglects such gradient terms |27| . 

We always use x for the transverse coordinate inside 
the reaction plane, with positive values in the direction of 
the impact parameter b, and y for the transverse coordi- 
nate perpendicular to b. (In momentum space y denotes 
the rapidity; which meaning is implied should follow from 
the context.) z points in beam direction. 



A. The equation of state 

The set ( [2.l| ) contains 5 equations for 6 unknown fields 
e,n,p,v x ,v y ,v z . To close the system one needs an equa- 
tion of state (EOS) which relates pressure, energy and 
baryon density. The EOS for strongly interacting mat- 
ter involves a phase transition from a hadron resonance 
gas (HG) phase to a color-deconfined quark-gluon plasma 
(QGP) phase. Like many others before (see, e.g., [28|,p9|) 
we accomplish this by separately constructing an EOS for 
a resonance gas (EOS H) and for the QGP phase (EOS I) 
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FIG. 1. The three equations of state discussed in the text, 
at vanishing net baryon density. 

and matching the two via the Maxwell construction, in- 
voking a bag constant B to describe the different vacuum 
energy in the two phases. EOS H is constructed from the 
contributions of all known hadron resonances of masses 
up to 2 GeV; their repulsive short-range interactions are 
parametrized via a mean-field potential V(n) = \Kn 2 
with K = 0.45 GeVfm 3 |j§. The QGP is described as an 
ideal gas of massless quarks and gluons (EOS I) inside a 
large bag with bag constant B. The latter is tuned to the 
desired phase transition temperature: B 1 / 4 = 230 MeV 
gives T c {n = 0) = 164 MeV at vanishing net baryon den- 
sity. EOS I is given by the simple equation p(e,n) = he 
or dp/de = independent of n. 

In order to investigate the influence of the phase tran- 
sition on the anisotropic transverse flow pattern, we stu- 
died separately the equations of state EOS H and EOS I 



as well as the combined equation of state EOS Q which 
includes the phase transition between the first two as ob- 
tained from the Maxwell construction. Comparisons to 
data are only performed for EOS Q. Figure [l] shows all 
three equations of state for vanishing net baryon density 
n = while Fig. || gives for EOS Q the pressure as a 
function of both e and n. 
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FIG. 2. The equation of state EOS Q with a quark-hadron 
phase transition. The pressure is shown as a function of en- 
ergy and net baryon density, e and n. For each value of n there 
exists a minimum energy density e m i n (n) with corresponding 
pressure p m in(e m i n , n); below e m i n the pressure is set to zero 
by hand. 



B. Reduction to 2+1 dimensions 

At high collision energies, relativistic kinematics and 
its influence on the particle production process implies 
longitudinal boost-invariance of the collision fireball near 
midrapidity p2fl . (Of course, near the target and projec- 
tile rapidities longitudinal boost-invariance is broken by 
the finite amount of total available energy.) As a re- 
sult, the longitudinal velocity field scales as v z = z/t, 
and it is convenient to use a coordinate system spanned 
by longitudinal proper time r = ty/l — v 2 z and space- 
time rapidity ri = h ln[(t+z)/(t— z)] instead of t and z 
(see Appendix |A|) . Longitudinal boost-invariance is then 
equivalent with 77-independcncc. 

Assuming the validity of this scaling ansatz near 
midrapidity, the longitudinal expansion of the fireball 
can be dealt with analytically, thereby reducing the nu- 
merical problem to the two transverse dimensions and 
time M. This greatly reduces the numerical effort. How- 
ever, by doing so one gives up the possibility to study the 
rapidity depen dence of the (anisotropic) transverse flow 
pattern PoPq ] as well as other interesting effects which 
occur at AGS and SPS energies, like the tilt of the longi- 
tudinal axis of the collision fireball away from the beam 
direction p6pC| Jl1]. For such studies a complete solution 
of the (3+l)-dimcnsional hydrodynamics |g^,^,^,^3) is 
required. We will here concentrate entirely on the midra- 
pidity region where the (2+l)-dimensional approach with 
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exact longitudinal boost-invariance is expected to yield 
reasonable results even at SPS energies. At higher ener- 
gies the model should become better and better. 

The implementation of longitudinal boost-invariance 
and transformation from (t, z) to (r, 77) is described in Ap- 
pendix The rewritten hydrodynamic equations read 



d T T TT + d x (v x T TT ) 

d T f™ + d x (v x f™j 
d T f r y + d x (v x f T yj 

drf + d x (v x j T ) 
where 



dy WUyT TT 
dy (Vyf™ 
dy (vyfTV- 

dy(VyJ T ) 



= ~P, 

= -d x p , 

= -d y p , 
= 0, 



(2.4) 



T» v = tT^ , p = 
Vi — Vi cosh r\ , Vi 



Tp, 

rji TT 



(e +p)7 2 — p 



(2.5) 
x,y)- 



We call Vi the transport velocities and Vi the energy flow 
velocities in the transverse directions. Since we work at 
midrapidity, 77 = 0, the transverse transport velocities 
agree with the corresponding fluid velocities in the cm. 
frame . 

In hydrodynamic problems phase transitions generi- 
cally lead to the formation of shock waves which compli- 
cate the n ume rical solution. To integrate the differential 
equations ( |2.4| ) we use the "Sharp and Smooth Transport 
Algorithm" (SHASTA |3l| ) which was shown to perform 
excellently even under difficult conditions |B3| . 



C. Initialization of the fields 

In this subsection we discuss the initial conditions for 
the solution of Eqs. (fO|). Strong interactions between 
the partons of the colliding nuclei lead to the deposition 
of a large fraction of the beam energy and the creation 
of many secondary particles in the reaction zone. The 
newly produced partons interact strongly with each other 
and, after only a few scatterings during a time interval 
To = 0(1 fm/c), the system is expected to reach a state of 
approximate local thermal equilibrium. Following [p|j32| 
(to which we refer for details) we take the energy depo- 
sition in the transverse plane to be proportional (by a 
factor K) to the number of collisions producing wounded 
nucleons: 



e{x,y;T ) = 
K\T A (x+ly) 



T B {x-±y) 



(2.6) 



, (1 oT B (x-\,y) y 



1-1- 



A 



Here Ta is the nuclear thickness function of the incoming 
nucleus A, 



+00 



T A (x,y) = / dzp A (x,y, z) 



(2.7) 



where the nuclear density pA is given by a Woods-Saxon 
profile, 



Pa{t) 



Pa 



1 + exp[(r - Ro)/(\ 



(2.8) 



and similarly for nucleus B. 

We further assume that the initial transverse density 
profile of net baryon number is proportional to the initial 
transverse energy density profile: 



n(x,y;T ) = Le(x,y;ro) ■ 



(2.9) 



For Pb-Pb collisions we use in ( |2.8| ) a nuclear radius Rq = 
6.5 fm and a surface thickness £ = 0.54 fm J34|. For U- 
U collisions we take Rq — 6.8 fm, with a deformation 
S = 0.27 ( H], Vol. 2, p. 133). This leads to a ratio 
Ri/R s = 1-29 between the long and short axes of this 
nucleus; in absolute terms Ri — 8.0 fm and R s = 6.2 fm 
p5[ . For the ground-state nuclear density we take p — 
0.17fm- 3 [51. 
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FIG. 3. Left: initial transverse energy density distribu- 
tion for a typical 158AGeV/c Pb+Pb collision at impact 
parameter b = 7 fm. Indicated are contours of constant ener- 
gy density between e = 7.0GeV/fm 3 (innermost contour) and 
e = 0.5GeV/fm 3 (outermost contour) in steps of Ae = 0.5 
GeV/fm 3 . The dashed lines represent the colliding nuclei 
before impact. Right: The same for a central 155AGeV/c 
side-on-side U+U collision - the innermost (outermost) con- 
tour corresponds to e = 8.0GeV/fm 3 (0.5 GeV/fm 3 ). 

Three parameters thus describe the initial conditions: 



• the maximum energy density eo in a central colli- 
sion (6 = 0)- this fixes the parameter K in (2.6) 
at the given beam energy; 



• the ratio L in ( |2.9|) between energy and baryon den- 
sity; 



the equilibration time tq. 
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In Sec. [IE we adjust the parameters by tuning the out- 
put of our calculations with EOS Q for central (b = 0) 
Pb+Pb collisions to experimental data (transverse mass 
spectra of negative hadrons and net protons at midrapid- 
ity pq] ) at 158 AGeV/c beam momentum. We use the 
same parameters K, L and To for U+U collisions at 155 
AGcV/c. 

In Fig. |^ we illustrate the initial conditions resulting 
from this tuning procedure. It shows contour plots of 
the energy density in the transverse plane at z = for 
Pb+Pb collisions with 6 — 7fm and central U+U colli- 
sions in the side-on-side configuration at the highest SPS 
beam momentum of 400 x (Z/A) GcV/c. Note that at 
fixed collision energy the central energy density for 6 = 
side-on-sidc U+U collisions is 8% lower than for 6 = 
Pb+Pb collisions, but about 14% higher than in Pb+Pb 
collisions at 6 = 7 fm which correspond to about the same 
initial spatial deformation. At similar deformation, the 
initial volume of the elliptic fireball formed in central 
side-on-side U+U collisions is almost twice that of the 
corresponding semi-central Pb+Pb collisions. 



D. Freeze-out and particle spectra 

As the matter expands and cools, the mean free path of 
the matter constituents grows, and the hydrodynamical 
description eventually breaks down. The system reaches 
the point of "kinetic freeze-out" after which the momen- 
tum spectra are no longer significantly affected by scat- 
tering among the particles. One should stop the hydro- 
dynamic solution when the average time between scatter- 
ings T sca tt = l/(va)n becomes comparable to the expan- 
sion time scale r cxp = 1/9 • u (inverse "Hubble constant" ) 
p9| [IH . [It was shown in that in relativistic heavy 
ion collisions freeze-out happens dynamically rather than 
geometrically, i.e. it is driven by the expansion of the 
fireball and not by its finite size.] Numerical calcula- 
tions |42|]43| l have shown that, since the particle density 
in the denominator of r scat t is a very steep function of 
T, this leads to freeze-out at nearly constant tempera- 
ture. For low net baryon freeze-out densities, as they 
arise in heavy ion collisions at and above SPS energies 
near midrapidity, this corresponds to almost constant en- 
ergy density. We here therefore impose freeze-out at a 
constant energy density edcc which is the most easily im- 
plemented condition in hydrodynamics. The value of edec 
(or, nearly equivalently, Td oc ) is another model parame- 
ter to be tuned to the data. 

After the freeze-out hypersurface £ of constant energy 
density edoc has been determined, the Td ec (a;)j chemical 
potentials Hi(x) and flow velocity field u^x) are eval- 
uated on this surface. To this end a tabulated version 
of EOS H is used for interpolation which (in addition 
to the pressure p) gives the intensive thermodynamical 
variables as functions of e and n. Each cell x on this 
freeze-out hypersurface contributes particles of species i 



(where i runs over all resonances included in EOS H) 
with a local equilibrium distribution 



fi(x,p) 



9i 



1 



(27T) 3 e lpMx)-Vi(x)]/T doc (x) ± I 



(2.10) 



<7i is the spin-isospin degeneracy factor for particle species 
i. The complete momentum spectrum is obtained by 
summing the corresponding particle flux currents across 
the 3-dimensional freeze-out hypersurface £ in space- 
time over all cells in £ (Cooper-Frye prescription |23|]): 



E 



,dNi 



dN, 



d 3 p dyp T dp T dip 



p-d 3 a(x)Mx,p). (2.11) 



This prescription is strictly correct only for freeze-out 
surfaces whose normal vector d 3 a(x) is everywhere time- 
like because otherwise some particles flow back into the 
4- volume inside £. A discussion of this issue which still 
awaits a fully consistent solution can be found in Q,(l5| . 

In the present paper we concentrate on flow patterns 
reflected in pion spectra. (Flow anisotropies f or p ions 
and protons at SPS energies were compared in Q.) A 
significant fraction of the measured pions arises from the 
decays of unstable resonances after freeze-out. These de- 
cays usually happen isotropically in the rest frame of 
the resonance and tend to smear out flow anisotropies, 
thereby reducing the anisotropic flow signals [p+|j2^ 1 . The 
fraction of pions from resonance decays depends strongly 
on the freeze-out temperature: their diluting effect on 
the elliptic flow V2, for example, is much stronger at 
Tdoc = 140 MeV fH than at T Acc = 120 MeV 0. All our 
calculations fully account for resonance decay contribu- 
tions including the complete relativistic decay kinematics 

M. 



E. Tuning the model 

Since the hydrodynamic approach cannot describe the 
initial thermalization stage directly after nuclear impact, 
the initial conditions for the hydrodynamic expansion 
stage cannot be predicted but must be obtained by fit- 
ting experimental data. However, once the initial con- 
ditions (in our case the parameters K, L, and r ) have 
been fixed in central collisions, the Glauber model ( |2.6| ) 
uniquely predicts their dependence on the impact param- 
eter. The validity of the hydrodynamic model can thus 
be tested by checking the impact parameter dependence 
of its predictions. In jll],[l7j we showed that, after be- 
ing tuned to central Pb+Pb collisions at 158 AGeV/c, 
the model successfully reproduces the measured pion and 
proton spectra near midrapidity up to impact parameters 
of 8-10 fm. This was better than expected. 

We here provide some details of the tuning procedure 
which were not previously reported in [p.l[ due to space 
limitations. In particular we show in Fig. |4| our fit to the 
midrapidity m T -spectra of negative hadrons (h~ ) and net 
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protons measured by the NA49 collaboration |3(|. The 
theoretical spectra are absolutely normalized. The corre- 
sponding fit parameters for the initial state are eo = 9.0 
GeV/fm 3 for the initial energy density in the cent er o f 
the fireball (corresponding to K — 2.04GeV/fm in (|2.6| ) 
and to an initial central temperature To = 258 MeV |37fl ), 
no = 1.1 fm~ 3 for the initial baryon density in th e fir e- 
ball center (corresponding to L = 0.122 GeV -1 in ( |2.9| )). 
and a starting time tq = 0.8fm/c for the hydrodynamic 
expansion (corresponding to T ■tq/%= 1.05). r controls 
the dilution of the matter via boost-invariant longitudi- 
nal expansion and thus the length of time available for 
the buildup of transverse flow before freeze-out; the latter 
affects the slope of the m T -spectra. The total time until 
freeze-out and the amount of transverse flow generated 
can also be changed by varying the initial energy density, 
but this also affects the normalization of the midrapidity 
spectra, eo and tq result from a suitable balance between 
these two effects, no is then essentially fixed by the mea- 
sured ratio between the proton and h~ spectra. 




FIG. 4. Particle spectra from central Pb+Pb collisions at 
158AGeV/c at midrapidity together with the hydrody- 
namical model predictions after tuning of the model parame- 
ters (solid lines). 

The different shapes of the proton and h~ spectra 
provide a handle to separate collective transverse flow 
(((v±)}) from thermal motion (Td cc ) at freeze-out. How- 
ever, it is known that a thermal model analysis of par- 
ticle spectra in general results in strong correlations be- 
tween these two parameters pl|,p2[ . Our best fit gives 
T doc wl20MeV (corresponding to e doc = 0.06 GeV/fm 3 ) 
and = 0.45 c, albeit with a significant uncertainty 

(somewhat lower Td ec with higher ((v±)) and vice versa 
cannot be excluded). This is in good agreement with 
other analyses of particle spectra Eq] and hydrodynamic 
simulations Ej^]; a combined analysis of spectral slopes 
and two-particle Bose-Einstein correlations |5(],|5l| tends 
to give somewhat larger transverse flow velocities coupled 
to lower freeze-out temperatures, but still inside the re- 
gion of uncertainty from the analysis of the single-particle 



spectra. 

This set of fit parameters, adjusted to SPS data, is our 
starting point for extrapolations towards non-central col- 
lisions and into different collision energy regimes. When 
studying the impact parameter dependence at fixed col- 
lision energy we leave all parameters unchanged. This 
may be unrealistic for very peripheral collisions where 
the midrapidity fireball is smaller and geometric freeze- 
out can cut the expansion short, leading to higher de- 
coupling temperatures. For the spectral slopes this is a 
second order effect since earlier freeze-out at higher T^ ec 
is partially compensated for by a smaller transverse flow 
velocity As we will see below (see Fig. ^), the 

elliptic flow anisotropy «2 builds up early in the colli- 
sion and, even at SPS energies, has almost reached its 
final value already several fm/c before decoupling; a pos- 
sible earlier decoupling in very peripheral collisions thus 
will not strongly affect either. We thus feel justified 
in leaving the model parameters (in particular the decou- 
pling temperature) unchanged when studying the impact 
parameter dependence. 

When investigating the excitation function of radial 
and elliptic flow we change K and tq. This is rationalized 
as follows: At higher energies we expect higher particle 
production per wounded nucleon; we cannot predict the 
beam energy dependence of secondary particle produc- 
tion, but we can parametrize it by changing K and plot- 
ting our results as a function of the finally observed mul- 
tiplicity density dN/dy. The beam energy dependence of 
dN/dy will eventually be provided by experiment, then 
allowing to present our results directly against ^/s. 
Higher initial particle production leads to higher particle 
and energy densities and thereby to accelerated thermal- 
ization. From relativistic kinematics and the uncertainty 
relation it follows that the production time of a secondary 
particle is inversely related to its energy |52j ; by dimen- 
sional analysis this suggests that the thermalization time 
To scales in inverse proportion to the initial temperature 
To: ToTo = const, or, equivalently, tq A' 1 / 4 = const. This 
is what we use in the present paper; in |Tl|] we instead left 
To constant. Within the range of collision energies stud- 
ied in Jll| the difference is negligible, but for the higher 
energies investigated here a reduction of tq ~ 1 /Tq causes 
a significant shrinkage of the horizontal axis on the exci- 
tation function in Fig. |TJ below. 

For energies above the SPS we leave the initial baryon 
density n(x, y; tq) unchanged. As a result, the ratio L 
of baryon to energy density drops, qualitatively consis- 
tent with the expectation of decreasing baryon stopping. 
Since already at the SPS the influence of the baryons 
on the EOS is minor, it doesn't really matter in which 
way L approaches zero as the collision energy goes to 
infinity. Note that we don't predict the normalization 
of the baryon spectra at other than SPS energies. Be- 
low SPS energies we leave L constant, lacking motiva- 
tion for a different choice. Once a better understanding 
of the beam energy dependence of the initial conditions 
becomes available, this can be easily improved. 
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III. TRANSVERSE FLOW PHENOMENOLOGY 

In this section we study generally the space-time evo- 
lution of the transverse flow pattern and how it is in- 
fluenced by a phase transition in the EOS. Since the 
finally observed particle spectra and their azimuthal 
anisotropies reflect the full space-time history of the fire- 
ball expansion, their proper interpretation requires an ac- 
curate understanding of the transverse fireball evolution. 
In jllj we showed that the softening of the EOS in the 
phase transition region leads at collision energies above 
the SPS to a reduction of the elliptic flow coefficient V2 
below the value expected from a hadron resonance gas. 
At even higher energies, however, one expects to enter a 
regime where the initial energy density is so far above the 
phase transition that nearly all of the expansion history 
happens inside the QGP phase. Since far above T c the 
EOS of a QGP (p = \e - B) is much harder than EOS H 
(which in the region relevant for us can be parametrized 
by p«0.15e), V2 should eventually rise again and ap- 
proach the value characteristic of EOS I which is 30-40% 
higher. In order to see whether this is true we have now 
studied collisions at very much higher energies, even far 
beyond the LHC. 



A. Semiperipheral Pb+Pb collisions 

In this subsection we investigate Pb+Pb collisions at 

an impact parameter of 7fm (left panel in Fig. ||). We 

begin by showing the evolution of the energy distribution 

and flow field in the transverse (x, y) plane for the cases 

with and without a phase transition. We do so for an 

initial central energy density in b — Pb+Pb collisions of 

e = 175 GeV/fm 3 (T = 510 MeV) at r = 0.38 fm/c. The 

resulting total pion multiplicity density with EOS Q of 

^r JL \ „ = 1070 at 6 = 7fm is at the upper end of the 
dy \ y =o *21 

range of predictions for RHIC energies |p3[ . This study 
was motivated by the work of Teaney and Shuryak who 
predicted under similar conditions an interesting phe- 
nomenon which they called "nutcracker flow" jjjj and 
which shows up only in the presence of a phase transi- 
tion. In Fig. |5| wc show the evolution for EOS I, i.e. a 
hard EOS without phase transition. One sees smooth ex- 
pansion and a continuous transition from an initial state 
of positive elliptic deformation (longer axis perpendicu- 
lar to the collision plane) to one with negative deforma- 
tion, caused by the developing in-plane elliptic flow. The 
thicker contours correspond (from the inside outward) to 
e = 1.6, 0.45, and 0.06GeV/fm 3 ; for the more realistic 
equation of state EOS Q the first two values limit the 
mixed phase while the latter indicates freeze-out. 

Figure || shows the analogous situation for EOS Q 
(which includes a phase transition) for identical initial 
conditions. Compared to Fig. ||one sees clear differences: 
the lack of a pressure gradient in the mixed phase inhibits 
its transverse expansion; the hadronic phase ouside the 



mixed phase expands quickly and freezes out, leaving a 
shell of mixed phase matter behind which inertially con- 
fines the QGP matter in the center. The matter with the 
softest EOS (smallest p/e) is concentrated around the 
QGP/mixcd interface (thick contour at 1.6 GeV/fm 3 ). 
When the QGP matter finally pushes the mixed phase 
shell apart (the "nutcracker phenomenon" discovered in 
|]l2| ), the energy density contours develop an interesting 
structure vaguely reminiscent of two separated half shells. 
Compared to Fig. ||, the elliptic flow clearly needs more 
time to push the matter from a state of positive to one of 
negative elliptic deformation. This is due to the inertia 
of the mixed phase shell which does not participate in 
the pushing. 




5 10-10 -5 5 10 
x(fm) x(fm) 
FIG. 5. Time evolution for EOS I of the transverse en- 
ergy density profile (indicated by constant energy density con- 
tours spaced by Ae = 150MeV/fm 3 ) and of the flow velocity 
field (indicated by arrows) for Pb+Pb collisions at impact pa- 
rameter 6 = 7.0 mi. The four panels show snapshots at times 
r— to = 3.2, 4.0, 5.6, and 8.0 fm/c. At these times the max- 
imal energy densities in the center are 5.63, 3.62, 1.31 and 
0.21 GeV/fm 3 , respectively. For further details see text. 



Figures || and ^| emphasize the spatial structure of the 
fireball at fixed time steps. Let us now study the time 
evolution in more detail. To this end we condense the 
information contained in the density and flow patterns 
into three time-dependent scalar quantities: 

(i) The "spatial ellipticity" 



«y 2 



((y 2 +x 2 )) 



(3.1) 



characterizes the spatial deformation of the fireball in 
the transverse plane. The angular brackets denote en- 
ergy density weighted spatial averages at a fixed time. 
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e x causes azimuthal anisotropies in the transverse pres- 
sure gradients which would eventually drive it to zero if 
the hydrodynamic evolution were not cut short by the 
freeze-out process. 




10-10 -5 5 10 

x(fm) x(fm) 
FIG. 6. Same as Fig. |5[ but for EOS Q which features a 
phase transition. The spacing between energy density con- 
tours is again 150 MeV/fm 3 , and the snapshots are taken at 
the same times. The corresponding maximum energy densi- 
ties are 5.97, 3.97, 1.67, and 0.55 GeV/fm 3 , respectively. See 
text for discussion. 



(ii) The momentum anisotropy 



(3.2) 



measures in an analogous way the anisotropy of the trans- 
verse momentum-space density. It is directly calculated 
from the spatial components of the energy momentum 
tensor but, as shown in [jOl , at freeze-out it is nearly 
equal to the p^,-weighted elliptic flow v 2 ^ for pions as 
calculated from their final momentum spectra ]5jJ |. Its 
time-dependence thus provides a picture of the dynami- 
cal buildup of the elliptic flow even at early times when 
the elliptic flow coefficient v 2 (which is calculated from 
hadronic momentum spectra, see Sec. IV) is not yet de- 
fined. For pions at freeze-out v 2 is given by 2 v 2 ~ 
V 2.,p% ~ e P |1- 

(in) The time-dependence of the average radial flow 
velocity 



l\]v 2 x +vl 
((7)) 



(3.3) 



the time-dependencies of {(v±)) and e p allows to answer 
the question to which stages of the expansion (i.e. to 
which domains of the EOS) each one is most sensitive. 
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characterizes the buildup of the overall transverse expan- 
sion which is modulated by the elliptic flow. Comparing 



FIG. 7. Time evolution of the spatial ellipticity e x , the mo- 
mentum anisotropy e p , and the radial flow {{v±)). The labels 
a, b, c and d denote systems with initial energy densities 
of 9, 25, 175 and 25000 GeV/fm 3 , respectively, expanding 
under the influence of EOS Q. Curves e show the limiting 
behaviour for EOS I as eo — > oo (see text). In the lower 
two panels the two vertical lines below each of the curves a-d 
limit the time interval during which the fireball center is in 
the mixed phase. In the upper panel the dots (crosses) indi- 
cate the time at which the center of the reaction zone passes 
from the QGP to the mixed phase (from the mixed to the HG 
phase). For curves a and b the stars indicate the freeze-out 
point; for curves c-e freeze-out happens outside the diagram. 



We now give a detailed discussion of Figs. [f]a-c which 
show (using EOS Q) the time evolution for the above 
three quantities for a sequence of collision energies, 
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parametrized by the initial central energy density in b = 
Pb+Pb collisions, e : e = 9, 25, 175, and 25000 GeV/fm 3 
(curves a through d in Figs. 0). With increasing eo the 
initial time tq was scaled down as described at the end 
of Sec. HE. The lowest of these eo-values corresponds to 
158AGeV Pb+Pb collisions at the SPS, while the high- 
est value is far beyond the reach of even the LHC. 

A calculation with EOS I is shown for comparison as 
curve e. Since EOS I (e = 3p) is completely scale in- 
variant, the time evolution of the dimensionless ratios 
( |3.ip , (3.2), and ( [3.3| ) is invariant und er a rescaling of 
eo as long as r is held fixed (see Eqs. (2.4)). Changing 
To ~ e breaks this scaling, but only weakly as we 
have checked. Curves e in Fig. show the time evolution 
for EOS I in the limit eo — > oo, to — * 0. Not shown is 
a calculation with EOS Q which was initialized with an 
extraordinarily high initial temperature of To « 20 GeV 
(eo = 25 x 10 6 GeV/fm 3 ); during the first 16 fm/c covered 
by Fig. [?] it fully coincides with curve e. In this case al- 
most all of the matter stays in the QGP phase during this 
time period whose EOS coincides with EOS I up to the 
(here negligible) bag constant. Therefore, as expected, 
the hydrodynamic evolution with EOS Q approaches at 
asymptotically high energies that with EOS I. 

Inspection of Fig. Q shows that the elliptic flow e p sat- 
urates at large times while the radial flow ((v±)) keeps 
rising forever, albeit at a decreasing rate. The driving 
force for radial flow, the radial pressure gradient between 
the matter in the fireball and the surrounding vacuum, 
never vanishes completely. The spatial ellipticity e x , 
on the other hand, which is responsible for azimuthal 
anisotropics in the transverse pressure gradients and thus 
drives the evolution of e p , passes through zero after some 
time. Afterwards the longer axis of the transverse fire- 
ball cross section no longer points perpendicular to the 
reaction plane, but into the reaction plane. A vanish- 
ing e x implies a vanishing growth rate for e p ; as e x turns 
negative, smaller oppositely directed anisotropies of the 
pressure gradients develop which can actually cause e p to 
decrease again. This can be seen in Fig. ^|b for large val- 
ues of eo where the sign of e x changes sufficiently early 
in the collision that pressures are still high enough to 
generate this effect. 

Qualitatively one hence can say that the final value 
of e p is established roughly at the point when e x passes 
through zero. For SPS energies this happens just before 
decoupling (implying that the fireball freezes out in a 
nearly circular configuration), but at high energies this 
occurs well before freeze-out. Generically the freeze-out 
value of e p (and thus V2) is sensitive to the EOS at signifi- 
cantly higher energy densities than the radial flow ((v±J) . 
The elliptic flow indeed measures the early pressure pig . 

On a more detailed level, the time evolution shows an 
interesting additional feature: In curves b and c the ellip- 
tic flow e p is seen to peak even before e x passes through 
zero. The origin of this phenomenon, which is rela ted to 
the phase transition, will be discussed in Sec. [II C below. 



Comparison of the lower two panels in Fig. [?] shows 
that the softening effect on the EOS of the phase transi- 
tion affects the buildup of ((v±)) and e p at similar times. 
However, the influence on e p is stronger since elliptic flow 
is a smaller effect (which feels only the anisotropies in the 
transverse pressure gradient, not its overall magnitude) 
and thus more fragile than radial flow. This results in 
a relatively larger sensitivity of elliptic flow to the phase 
transition. 



B. 



Central U+U collisions in the side-on-side 
configuration 



As discussed in Sec. II C, central U+U collisions in 



the side-on-side configuration provide 14% higher energy 
density over nearly twice the volume at the same ini- 
tial spatial deformation as Pb+Pb collisions at b = 7 fm. 
This leads to a longer lifetime for non-zero spatial ellip- 
ticity e x , the driving force for elliptic flow, and also for 
the whole fireball until freeze-out. Hence the system has 
more time for thermalization, favoring the applicability 
of our hydrodynamic method. For this reason we decided 
to perform quantitative calculations for this system and 
make predictions for experiments with uranium beams at 
RHIC and LHC. 
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FIG. 8. Same as Figs. | and | (eo = 175 GeV/fm 3 at 
to = 0.38 fm/c, EOS Q), but for central side-on-side U+U col- 
lisions. The spacing between energy density contours is again 
150 MeV/fm 3 , and the snapshots are taken at the same times. 
The corresponding maximum energy densities are 8.71, 6.06, 
3.27, and 1.47 GeV/fm 3 , respectively. See text for discussion. 



We first look once more at the space-time evolution of 
the transverse energy density and flow profiles, shown 
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in Fig. The initialization corresponds to the same 
collision energy as in Fig. |6| (eo = 175 GeV/fm 3 ) but, 
since we now consider central (b = 0) collisions, the ini- 
tial energy density in the center of the deformed collision 
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FIG. 9. Same as Fig. fj], but now comparing central U+U 
(solid) to semiperipheral (6=7 fm) Pb+Pb collisions (dashed) 
at two selected beam energies. The curves labelled "SPS" cor- 
respond to e = 9 GeV/fm 3 (8.3 GeV/fm 3 ) for central Pb+Pb 
(side-on-side U+U) collisions, those labelled "RHIC" have 
eo = 175 GeV/fm 3 in both cases. 



region is higher than in the semiperipheral Pb+Pb col- 
lisions of Fig. [| As seen in Fig. [|, the whole time evo- 
lution is slower for central U+U than for semiperiphe- 
ral Pb+Pb collisions, due to the larger system size: At 
t — t o = 3.2 fm/c (the first shown snapshot) the central en- 
ergy density is 50% higher, and at t— tq = 8 fm/c (the last 
snapshot) it is even by a factor 3 larger than in b = 7 fm 



Pb+Pb collisions at the same beam energy. Freeze-out 
occurs nearly 30% later in central U+U than in semipe- 
ripheral Pb+Pb collisions (see Fig. |]). 

We note with surprise that the "nutcracker" phenome- 
non Il2] is conspicuously missing in the U+U collisions. 
We could not find it at lower and higher collision energies 
either. The origin of this difference between central U+U 
and peripheral Pb+Pb collisions will be discussed in the 
following subsection. 

In Fig. H we compare the time evolutions of the 
three characteristic quantities e x , e p , and ((v±)) in cen- 
tral U+U and semiperipheral Pb+Pb collisions, at SPS 
(e = 9GcV/fm 3 ) and RHIC (e = 175 GeV/fm 3 ) ener- 
gies. We note that at freeze-out (Td ec = 120 MeV) both 
systems give nearly the same radial and elliptic flow, in 
spite of the different time evolution: in the large system 
both flow types develop more slowly, but over a longer 
time. This does not take into account that the flow gra- 
dients are smaller in the larger system, leading to later 
freeze-out at a lower temperature p9| ]. This would not 
change the elliptic flow since e p has already saturated 
(actually, it would lead to a very slight decrease of e p , 
see Fig. ^|b). The radial flow ((v±)) would, however, be 
somewhat larger. Since we enforced freeze-out at the 
same value T^ ec , we don't see this. 



C. What makes the nut crack? 

In this subsection we analyze two questions which so 
far remained open: (1) Why does the "nutcracker" phe- 
nomenon arise in semiperipheral Pb+Pb collisions, but 
not in central U+U collisions, in spite of their identical 
initial deformation? (2) What is the origin of the de- 
crease of e p (r) before e x passes through zero which is 
observed in Fig. and curves b and c of Fig. fj]b? 

To answer them requires a more detailed look at 
the time evolution of the transverse pressure gradients 
(cause) and transverse flow profiles (effect). In Figs, [h] 
and [Til we show a series of six snapshots each for semipe- 
ripheral Pb+Pb and central U+U collisions, plotting the 
pressure and flow velocity profiles along the x and y axis, 
respectively. The crucial difference between the two col- 
lision systems is that in the semiperipheral Pb+Pb col- 
lisions the initial fireball contains a roughly 0.5 fm thick 
layer of mixed phase matter with vanishing transverse 
flow velocity; for central U+U collisions the initial ener- 
gy density drops to zero so steeply that the mixed phase 
layer is initially practically absent. 

As the matter begins to expand and dilute, a mixed 
phase layer begins to develop also in the U+U collisions; 
however, due to the buildup of transverse flow in the ex- 
panding matter, it is automatically created with a non- 
vanishing transverse flow velocity. Thus, even without 
pressure gradients inside the mixed phase which could ac- 
celerate it, the mixed phase matter flows in the transverse 
directions, with velocities exceeding those of the enclosed 
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QGP matter (see Fig.[Tl|). The resulting transverse flow 
profiles are monotonous functions of x and y, with a self- 
similar (linear "scaling" ) pattern inside the mixed phase 
exactly as given by the analytic solution recently found 
by Biro [ p9| . The monotony of the transverse flow pro- 
files is related one-to-one to the absence of the nutcracker 
phenomenon. 
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FIG. 10. Transverse pressure (solid) and velocity (dashed) 
profile, in x (thick) and y (thin) directions, for Pb+Pb col- 
lisions at b = 7 fm. The 6 panels show snapshots at the indi- 
cated times. The region of nearly constant pressure is in the 
mixed phase. The velocity profiles (dashed) are cut off at the 
freeze-out point. Initial conditions as in Fig. pi 



In the semiperipheral Pb+Pb collisions, on the other 
hand, the initially present mixed phase layer is at rest 
and, due to the lack of pressure gradients, cannot accel- 
erate itself in the transverse direction. As the transverse 
pressure gradients in the enclosed QGP matter begin to 
accelerate the QGP matter, the latter "slams" into the 
motionless mixed phase. This is clearly seen in the first 
four panels of Fig. ^ which show a strong radial increase 
of the transverse flow velocities inside the QGP phase, 
followed by a dramatic drop inside the mixed phase and 
a second rise in the HG matter near the edge. Inside 
the mixed phase the radial velocity profile is thus com- 
pletely different from the selfsimilar scaling pattern seen 
in Fig.|ll|. As time proceeds, this anomalous structure 
in the Pb+Pb collisions weakens, and the velocity profile 
begins to approach a scaling form inside the mixed phase; 
scaling violations survive longest near the outer edge of 



the mixed phase layer. In the y direction they disappear 
slightly earlier than in the shorter x direction; this is the 
origin of the "nutcracker phenomenon" . 
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FIG. 11. Same as Fig. [jij, but for central U+U collisions in 
the side-on-side configuration. Initial conditions as in Fig. ^. 

Now we can also understand the decrease of e p even be- 
fore e x passes through zero: Figs.]?], [lC] and [Tl| show that 
this happens while most of the fireball is in the mixed 
phase. (Actually, e p begins to decrease while there is 
still a small QGP core in the center.) During this stage 
the matter expands essentially without transverse accel- 
eration, featuring a nearly selfsimilar transverse flow pat- 
tern. While it lasts, the selfsimilar flow dilutes the ear- 
lier developed momentum anisotropy e p . This feature is 
therefore also directly related to the phase transition. 



IV. EXPERIMENTAL PREDICTIONS 

While the time-evolution of e x , e p and ((v±)) is interest- 
ing and helpful for an understanding of the relevant phys- 
ical mechanisms, only the final values at freeze-out are 
observable (through the momentum spectra and, in the 
case of e x , possibly indirectly via two-particle momentum 
correlations). The flow observables thus represent time- 
integrals over the expansion history and EOS, and their 
measurement in a single collision system at fixed beam 
energy provides very little information. Using flow signa- 
tures as indicators for properties of the equation of state 
for strongly interacting matter requires their measure- 
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ment over a wide range of external control parameters, 
such as impact parameter, size of the colliding nuclei, 
and beam energy. As discussed in the preceding section, 
a time-differential measurement is to some extent possi- 
ble by comparing the radial and elliptic flow as functions 
of these parameters. 

Flow anisotropies reflect themselves as non-vanishing 
higher order Fourier coefficients in a Fourier expansion of 
the azimuthal dependence of the measured single-particle 
spectra around the beam direction |p5[ : 

f 77 dip cosinip) 

I \ J ~ Tr ay dtp , tAi\ 

Vn{y) = fS— — dN ' n = l,2,.... (4.1) 



dy dip 



Since most experiments have limited p T acceptance, one 
studies these coefficients also as functions of the trans- 
verse momentum: 

f n dip cos(nip) — 2d™ — _ 
v n (y,P T ) = — I , *L*p&z*r . (4.2) 



dN 
I pt dpT dip 



The p 2 -weighted anisotropic flow coefficients are defined 
by 



J^dip cos(^) $ v \d v \ 1 -^ 



Vn,p 2 T (y) = 71T-, f 

i-» dl P J 



dy dp-*, dip 



FT rT dy dp? v d(p 



(4.3) 



In symmetric collision systems (which are the only ones 
we consider here) the odd order coefficients V\ , 1)3, . . . van- 
ish at midrapidity y = by symmetry. We here concen- 
trate on the second harmonic coefficient which is conven- 
tionally called "elliptic flow" . The Vi are only defined at 
freeze-out but we already discussed how v 2 and v 2p -^ can 
be related to e p which is known also before freeze-out. 

e x and e p are functions of time; in the present sec- 
tion, however, we only need the initial spatial deforma- 
tion e x (to) and the final momentum-space deformation 
e p (jt). For simplicity we will quote them as e x and e p , 
respectively, without the time arguments. 



A. pT-dependence of elliptic flow 

Since most experiments have a limited acceptance in 
transverse momentum, the measured elliptic flow signal 
must be corrected for the p T -acceptance. In Fig. |l2| we 
show the p T -dependence of v 2 for pions and protons for 
scmipcripheral Pb+Pb and central U+U collisions. In 
spite of their different masses, the predicted V2(p T ) is 
rather similar for the two particle species |Tcfl . At low 
p T , the heavier protons show even a little less elliptic flow 
than the pions. To the extent that hydrodynamics is ap- 
plicable, the larger (v 2 ) for protons than pions measured 
by NA49 [|f7| is thus predominantly due to the differ- 
ent p T -windows for the two particle species (the proton 
elliptic flow was measured at higher p T JTt|] ) . 



According to general arguments [Q, v 2 must vanish 
with zero slope as p T — > 0. We checked that this is true. 
Fig. [l^ shows, however, that for pions the turnover from 
a roughly linear behaviour at large p T to zero slope as 
p T — * occurs at very small p T -values, p T < 0.1 GeV/c; 
for protons the corresponding scale is somewhat larger. 
We have no quantitative analytic understanding of this 
momentum scale but note that qualitatively similar be- 
haviour was found in [ p7[ using the kinetic UrQMD 
model. 



Pb+Pb, b=7 fm 



U+U, b=0 fm, side-on-side 




1 2 3 "0 1 2 3 

p T (GeV) p T (GeV) 

FIG. 12. p T -dependence of the elliptic flow coefficient 112 for 
pions (solid) and protons (dashed), for 158 AGeV/c Pb+Pb 
collisions at & = 7fm (left panel) and 155 ^4GeV/c U+U col- 
lisions at b = in the side-on-side configuration (right panel) . 



B. Impact parameter dependence of elliptic flow 

As one changes the impact parameter, the initial spa- 
tial deformation e x of the transverse cross section through 
the reaction zone varies as shown in Fig. 3 of Ref. [0. 
The stronger the initial ellipticity, the stronger is the 
hydrodynamic response to it, i.e. the larger are v 2 or 
e p at freeze-out. Ollitrault || showed that for an EOS 
with a constant velocity of sound, de/dp — const., the 
ratio e p /e x or, equivalently, v 2 jt x is independent of the 
impact parameter [Q. (Ollitrault Js] used the variable 
«2 !p 2 which is closely related to e p |54|. For pions v 2 and 
e p are related by a factor 2 (ll).) This scaling is bro- 
ken only for very peripheral collisions which freeze out 
before the elliptic flow builds up and saturates; thus in 
hydrodynamics v 2 / e x is constant over most of the impact 
parameter range. 

A phase transition is characterized by a strong drop 
of the sound velocity in the critical region (for a first or- 
der phase transition the sound velocity vanishes in the 
mixed phase). It is therefore interesting to reinvestigate 
the impact parameter dependence of v 2 /e x in the pres- 
ence of a phase transition. The impact parameter not 
only controls the initial spatial ellipticity of the fireball, 
but also (with less variation) its initial energy density. At 
a given beam energy, it is therefore possible to probe the 
EOS over a range of energy densities by varying the im- 
pact parameter. For a beam energy, at which in central 
collisions the initial energy density is not too far above 
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the phase transition, it may thus be possible to study 
the effect of the reduced sound velocity near the phase 
transition on the elliptic flow by changing the impact pa- 
rameter. Weak structures in Fig. 9 of Ref. S first indi- 
cated that the quark-hadron phase transition might thus 
become visible. Our analysis improves on that analysis 
by including resonance decays which tend to dilute the 
elliptic flow signature |TT| . 

In Fig. |l3] we study the impact parameter depen- 
dence of V2/e x in Pb+Pb collisions for three different 
initial central energy densities: e = 25 GeV/fm 3 (corre- 
sponding to a low energy run at RHIC), eo — 9 GeV/fm 3 
(corresponding to collisions at the highest SPS energy 
of 158 A GeV), and eo = 4.5 GeV/fm 3 (corresponding to 
lower SPS energies around 40 A GeV) . The calculated 
total pion multiplicity densities at b — and midrapid- 



ity are ^ | (6 = 0) = 859, 460, and 317, respectively. 
For later comparison with Fig. 14 we also quote the cor- 



responding rapidity densities for semiperipheral Pb+Pb 
collisions: ^| s=0 (6 = 7fm) =415, 220, and 148, respec- 
tively. 
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FIG. 13. The ratio of the elliptic flow coefficient V2 and the 
initial spatial ellipticity e x as a function of impact parameter 
b for Pb+Pb collisions. Results for three values of the initial 
central energy density at 6 = (eo = 4.5, 9.0 and 25 GeV/fm 3 ) 
are shown. Note the suppressed zero on the vertical axis. 



Fig. 13 shows that, at moderate impact parameters, 
the largest elliptic flow is generated at the lowest of these 
three beam energies. At very large impact parameters 
(where hydrodynamics becomes doubtful) the generated 
elliptic flow naturally drops to zero, since the overlap re- 
gion and its initial energy density are then too small and 
the matter freezes out before flow can develop. What is 
interesting, however, is that at higher beam energies the 
elliptic flow starts out lower than at eo = 4.5 GeV/fm 3 , 
but then Vijt x rises with increasing b. In fact, for 
eo = 9 GeV/fm 3 this ratio reaches at b = 11 fm nearly the 
same value as for central collisions at cq — 4.5 GeV/fm 3 . 

The decrease with rising beam energy of vijt x at mod- 
erate impact parameters was found |ll| to result from the 



softening of the EOS in the phase transition region. The 
soft matter near the transition point inhibits the buildup 
of elliptic flow. Going at fixed beam energy to larger im- 
pact parameters is like going at fixed impact parameter 
to lower beam energies: in both cases the initial energy 
density in the collision zone is reduced, and eventually the 
matter is dominated again by the relatively hard hadron 
gas. When read from right to left, the curves in Fig. [l3| 
can thus be viewed as different projections of the exci- 
tation function of elliptic flow which will be discussed 
below. We emphasize in particular the rise of V2/t x to- 
wards larger impact parameters at the high SPS and the 
low RHIC energy: without a phase transition this fea- 
ture would be absent. Unfortunately, these variations are 
small (at the level of a few percent), and very accurate 
measurements are required to identify them. 

Preliminary data from 158 A GeV Pb+Pb collisions 
p7[ show a monotonous decrease of 1)2 /e x with increas- 
ing impact parameter, instead of the nearly constant be- 
haviour predicted by hydrodynamics (see Fig. [l^). For 
b ^ 0, however, the data seem to approach the hydro- 
dynamic prediction. It is possible that semiperipheral 
Pb+Pb collisions do not equilibrate quickly enough to 
permit the elliptic flow to fully reach the hydrodynamic 
limit. Indeed, kinetic simulations with the RQMD code 
H[l7||24), where the collision centrality is coupled to the 
degree of local thermalization, are able to qualitatively 
explain the observed decrease of Vijt x with increasing 
impact parameter: more peripheral collisions lead to less 
equilibration and hence to a weaker elliptic flow response 
to the initial spatial ellipticity. When RQMD is modified 
to simulate an EOS with a quark-hadron phase transi- 
tion ||, the same generic decrease is superimposed on 
the rise of V2/e x at large b shown here (middle curve in 
Fig. 13); this results in a decrease of {v2/e x )(b) which is 
first steep, then flattens, then finally steepens again ||. 

It is evident that a proper understanding of the inter- 
esting features in the impact parameter dependence of 
V2/€x predicted in || for Pb+Pb collisions require the 
separation of pre-equilibrium effects from those induced 
by the softening of the EOS near the phase-transition. 
A collision system which is large enough to ensure suffi- 
ciently rapid thermalization for hydrodynamics to apply 
would make life much easier. We therefore suggest to 
study elliptic flow in side-on-side U+U collisions at zero 
impact parameter and search for the hydrodynamically 
predicted phase transition signatures in the beam energy 
dependence of elliptic flow. 



C. Beam energy dependence of elliptic flow 

The time-dependence of the flow patterns discussed in 
Sec. Ill reflects itself also in the beam energy dependence 
of elliptic flow. We already noted in |ll|] that the phase 
transition causes an non-monotonic excitation function 
for the elliptic flow coefficient V2- as the collision energy 
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is increased, vi first rises (at low energies the fireball 
freezes out before the elliptic flow can saturate) but then 
decreases again as the initial energy density rises above 
the QGP threshold. Wc now understand that this de- 
crease is intimately connected to the diluting effects of 
the selfsimilar fireball expansion in the mixed phase, even 
before the spatial defo rmat ion e x passes through zero (see 
the discussion in Sec. [II C .) Without a phase transition 
(EOS H) this does not happen (see dash-dotted lines in 
Fig. |l4|) ; the slight decrease of V2 with EOS H at asymp- 
totically high energies has a different origin, namely a 
reduction of e p by the opposite sign of the spatial fireball 
anisotropy after e x has passed through zero. 
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The comparison of semiperipheral Pb+Pb collisions 
with central U+U collisions in the upper and lower pan- 
els of Fig.|lJ shows that this non-monotonic behaviour of 
the excitation function for V2 is not sensitive to the exi- 
stence of the "nutcracker phenomenon" : the decrease of 
V2 below its maximum in the SPS regime is only slightly 
weaker in the U+U case than for Pb+Pb, although only 
the latter features a "cracking nut" . Since elliptic flow is 
a fragile phenomenon which is quite sensitive to incom- 
plete thermalization, we believe that the most promising 
route towards experimental verification of the phase tran- 
sition signature suggested here is to study the excitation 
function of V2 in largest available deformed collision sys- 
tem, namely central side-on-side U+U collisions. 

In Ref. ill we missed the fact that at asymptotically 
high energies the elliptic flow coefficient V2 must approach 
the larger value corresponding to the stiffer QGP equa- 
tion of state EOS I. We calculated in jllj the excitation 
function for b = 7 fm Pb+Pb collisions only up to mul- 
tiplicity densities ^f- 1 = 500 and concluded prema- 
turely that t>2 saturates at high collision energies at a 
value below the value corresponding to EOS H. Fig.|l4| 
extends the excitation functions for both Pb+Pb and 
U+U collisions to LHC energies and demonstrates that 
i>2 begins to rise again, eventually approaching the EOS I 
limit. The dip, which indicates the presence of the phase 
transition, thus only covers the energy range between 
SPS and RHIC. Note that in the same energy region 
also the radial flow {(v±J) (dashed lines in Fig.|l4|) is pre- 
dicted to grow more slowly with \fs than at lower and 
higher beam energies where the expansion is predomi- 
nantly driven by pure HG or pure QGP matter. 



D. Elliptic flow as an estimator for the 
thermalization time scale 
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FIG. 14. Excitation function of the elliptic flow coefficient 
«2 (left vertical axis) and the radial flow {{v±))/c (right verti- 
cal axis) , for Pb+Pb collisions at b = 7 fm (upper panel) and 
side-on-side U+U collisions at b — (lower panel) . The hori- 
zontal axis gives the total pion multiplicity density at midra- 
pidity, _ j as a measure for the collision energy. Hori- 

zontal arrows indicate the regions covered by SPS, RHIC, and 
LHC. In the lower panel LHC would start around 5000. 



Throughout this paper we have assumed early thermal- 
ization followed by hydrodynamic expansion. For a given 
initial deformation of the collision zone in the transverse 
plane (which can be calculated from geometry once the 
impact parameter is known, for example by a measure- 
ment of the number of spectator nucleons), this guaran- 
tees the maximum possible momentum-space response in 
the form of elliptic flow. Any delay in the thermalization 
process will lead to a reduction of the elliptic flow: even 
without secondary collisions the spatial deformation of 
the region occupied by the produced particles decreases 
by free-streaming, and if thermalization effectively sets in 
later, the resulting anisotropies in the pressure gradients 
will be smaller, leading to less elliptic flow. 

We can use the above demonstrated fact that, up to 
variations of the order of 20%, the hydrodynamic re- 
sponse V2 to the elliptic spatial deformation at thermal- 
ization is essentially constant: v^ y °/e x ~ const, ss 0.25. 
This allows to interpret the measured V2 in terms of an 
effective initial spatial deformation at the point of ther- 
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malization, i.e. at the beginning of the hydrodynamic 
evolution. It is clearly not a good approximation to ideal- 
ize the initial kinetic equilibration stage of the collision by 
a stage of collisionless free-streaming followed by hydro- 
dynamic expansion, thereby assuming a sudden, but de- 
layed transition from a non-equilibrium initial state to a 
fully thcrmalizcd fluid. Still, this simple-minded picture 
can be used to obtain a rough first order-of-magnitude 
guess of the thermalization time scale, based on a mea- 
surement of v 2 . 

To this end we note that under free-streaming the 
phase-space distribution evolves as 



f(r,p,t) = f [r- j^(t-t ),p,t 



(4.4) 



Using a Gaussian parametrization for the initial phase- 
space distribution of produced secondary particles, 



f(r,p,t ) = exp 



2Rl 



y 

2RI 



P X +Py 

2A 2 



(4.5) 



one easily finds 



£x0) = 



/ cPrr 2 cos(2(y9 r ) J d 3 pf(r,p,t) 



J cPrr 2 f d 3 pf(r,p,t) 



R, 



R l 



R 2 x + R 2 y + 2(cAt) 2 ' 



(4.6) 



where At — t—tg is the time delay between particle for- 
mation and thermalization. Assuming that e x (to+At) 
can be obtained from the measured v 2 by dividing by 
w 0.25, we can extract At by rewriting ([Oj) as 



e x (t + At) 
e x (to) 



1 



(cAtf 



R 2 {\ + 5 2 ) 



(4.7) 



where 8 parametrizes the initial deformation via R x — 
R(1S), R y = R(l+S) such that e x (* ) = 2S/(l+5 2 ). 

Inserting appropriate values for R and 5 one finds that 
for Pb+Pb collisions at b = 7 fm a dilution by 50% of 
the elliptic flow signal by initial free-streaming requires a 
time-delay of order 3.5 fm/c until thermalization sets in; 
for central U+U collisions in the side-on-side configura- 
tion At 5 fm/c of approximate free-streaming would 
be required to dilute the elliptic flow signal by 50%. This 
(admittedly rough) exercise demonstrates two points: (i) 
U+U collisions provide the better chance to observe the 
full hydrodynamic elliptic flow signal, and (ii) the ob- 
servation of less elliptic flow than hydrodynamically ex- 
pected can be used to obtain a rough estimate of the 
thermalization time scale in the initial collision stage. 



V. SUMMARY 

On the basis of hydrodynamic simulations we analyzed 
the sensitivity of radial and elliptic transverse flow at 



midrapidity to the quark-hadron phase transition. We 
modelled this phase transition as a strongly first or- 
der phase transition with a latent heat of about 1.15 
GeV/fm 3 . It manifests itself dynamically as an expan- 
ding shell of mixed phase matter inside which all pressure 
gradients and thus all hydrodynamic acceleration forces 
vanish. Compared to the situation of a pure HG or a pure 
QGP phase this leads to a reduction of both radial and el- 
liptic flow. Elliptic flow, as the more fragile phenomenon 
which is generated only by azimuthal anisotropics in the 
pressure gradients, shows a larger sensitivity to the phase 
transition than radial flow. Also, since we showed that 
it saturates well before freeze-out, it more directly re- 
flects the EOS during the early and dense stages of the 
expansion. 

As a tell-tale signature for the phase transition we pre- 
dict a non-monotonic excitation function for the elliptic 
flow coefficient v 2 as shown in Fig. [l4| In the present 
paper we explored in great detail the origin of the dip in 
v 2l which we predict to occur in the energy region be- 
tween the SPS and RHIC, by performing a careful ana- 
lysis of the space-time evolution of the anisotropic trans- 
verse flow pattern for a variety of collision energies. As 
the dynamical origin of the phase transition signature 
in v 2 we identified the existence of a large subvolume of 
mixed phase matter which undergoes nearly selfsimilar, 
acceleration-free expansion while it lasts. In addition to 
the v 2 excitation function it leaves traces in the impact 
parameter dependence of the response v 2 / e x of the ellip- 
tic flow to the initial spatial deformation of the collision 
zone, and in the (not directly measurable) time evolution 
of the flow anisotropy e p . 

When colliding spherical nuclei with each other, the 
measurement of elliptic flow requires selecting collisions 
at rather large impact parameters (b > 5 fm) in order to 
achieve a sufficiently large initial spatial deformation of 
the nuclear overlap region (reaction zone). Correspond- 
ingly the overall size of the elliptically deformed, expand- 
ing fireball is small, and one may doubt the applicability 
of our hydrodynamic approach. We here point out that 
central U+U collisions in the side-on-side configuration 
provide nearly twice larger collision volumes at similar 
deformation as Pb+Pb collisions at b = 7 fm and should 
thus exhibit hydrodynamic behaviour much more clearly. 

We therefore carefully compared central side-on-side 
U+U collisions with semiperipheral Pb+Pb collisions at 
all collision energies. We showed that the phase transi- 
tion signature in the v 2 excitation function manifests it- 
self similarly in both collision systems. The U+U system 
should thus be preferred for its presumed better hydrody- 
namical behaviour and for the larger particle multiplici- 
ties which improve the statistics of elliptic flow measure- 
ments. The phase transition signal appears to be slightly 
stronger in the smaller Pb+Pb system; we were able to 
trace this to the "nutcracker phenomenon" of Shuryak 
and Teaney ]l2[| which, unfortunately, only occurs in the 
Pb+Pb system. In trying to understand the fragility of 
"nutcracker flow" we found that it crucially relies on the 
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existence of a rather thick shell of mixed phase matter 
at rest in the initial state of fireball expansion, which 
surrounds a significant core of QGP. In response to in- 
ternal pressure gradients the QGP core starts to expand 
and "slams" into the surrounding shell of mixed phase at 
rest. This cannot happen in central U+U collisions since 
there the initial transverse energy density profile drops to 
zero so steeply that no visible mixed phase shell forms. 

We thus conclude that the interesting "nutcracker 
flow" phenomenon constitutes a very fragile variant of 
anisotropic flow which is not generated in central U+U 
collisions. If the fireballs formed in semiperiphcral 
Pb+Pb collision should turn out to be too small to 
achieve sufficient local thermalization for hydrodynamics 
to work, it may be unmeasurable. Fortunately, the ellip- 
tic flow signature for the phase transition is more robust 
and does not require the actual "cracking of the nut" ; it 
should be clearly visible in central U+U collisions. 

This raises the question how to experimentally select 
the side-on-side collision geometry. By requiring zero 
spectators one can trigger on configurations in which 
the colliding nuclei overlap completely in the transverse 
plane. This still allows for arbitrary, but (up to a sign) 
equal angles (#i = ±62) between the beam direction and 
the long axes of the two deformed nuclei. The interesting 
side-on-side configuration corresponds to 0\ — O2 = 90°. 
Since this configuration has the largest initial spatial de- 
formation in the transverse plane, it generates the largest 
elliptic flow V2', therefore, Shuryak [Q suggested a cut 
on large v% to select the side-on-side collision geometry. 
Unfortunately, the event-by-event fluctuations of V2 are 
so large that this off-line trigger is not expected to be 
very efficient p0| ; furthermore, it would introduce an in- 
convenient trigger bias into our suggested investigation 
of the dependence of V2 on various control parameters. 

We have not been able to come up with a more ef- 
ficient selection criterium. We checked that with ini- 



tial conditions calculated according to (2.6), the pro- 
duced charged particle multiplicity densities at midra- 
pidity vary by less than 5% between tip-on-tip and side- 
on-side collisions (with side-on-side collisions producing 
more particles, with slightly smaller (px) at freeze-out). 
Again this difference is well below the expected level of 
event-by-event fluctuations. I ts sm allness is explained by 
the fact that with the Ansatz ( |2.6| ) the amount of entropy 
dS/dy stopped at midrapidity is essentially independent 
of the orientation 9\ = ±82 (for 0° it is 1.3% larger than 
for 90°), and boost-invariant longitudinal expansion con- 
serves dS/ dy. At higher collision energies minijet produc- 
tion may overtake the soft particle production processes 
implicitly assumed in (2.6); instea d of scaling with the 
number of wounded nucleons as in (2.6), minijet produc- 
tion scales with the number of nuclcon-nucleon collisions, 
involving the product rather than the sum of the nuclear 
thickness functions appearing in (2.6). In this case tip- 
on-tip collisions are expected to generate considerably 
more entropy in the transverse plane at midrapidity than 
side-on-side collisions, and one could trigger on the lat- 



ter by selecting for zero spectators combined with low 
dN/dy(y=0). - In the absence of an efficient trigger 
for side-on-side U+U collisions at present-day collision 
energies one will be forced to compare with data which 
are averaged over all orientations 6± = ±62- The compu- 
tation of an orientation-averaged excitation function for 
V2 is, however, numerically expensive; we therefore post- 
pone it until experiments involving U+U are approved. 

Our prediction of a dip in the excitation function of i»2 
at midrapidity is directly related to the one by Rischke et 
al. Q] of a dip in the excitation function for directed flow 
at forward and backward rapidities: both rely on the soft- 
ening of the EOS near the phase transition which results 
in reduced hydrodynamic pressure gradients. We point 
out, however, that, as the collision energy increases, the 
time interval during which directed flow is generated (the 
nuclear transition time) becomes shorter and shorter, and 
the prospects for sufficiently fast local thermalization to 
validate hydrodynamic concepts thus become worse and 
worse. The opposite is true for elliptic flow: Figs. [?] and || 
show that the time interval over which elliptic flow builds 
up approaches at high collision energies a finite limit 
of about 7fm/c for semiperipharal Pb+Pb and about 
12fm/c for central U+U collisions. The density of pro- 
duced particles, on the other hand, continues to increase, 
leading to shorter and shorter thermalization times. The 
hydrodynamic description of elliptic flow buildup should 
thus become better with increasing collision energy. 

We finally comment on the sensitivity of the proposed 
phase transition signature to our simple modelling of the 
phase transition: we used a Maxwell construction be- 
tween the HG and QGP equations of state, leading to 
a strong first order phase transition with large latent 
heat. We don't believe that smoothing the phase tran- 
sition to a rapid crossover will qualitatively alter our re- 
sults: the only major change will be a replacement of 
the acceleration-free mixed phase by a transition region 
with non-zero, but nevertheless small pressure gradients. 
However, since elliptic flow signals are generically weak 
and the predicted effects from the phase transition are 
at a level of only about 10% of this signal, further hy- 
drodynamic simulations using a more realistic modelling 
of the EOS may be required for a reliable quantitative 
assessment of the expected experimental signal. 
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APPENDIX A: IMPLEMENTATION OF BOOST 
INVARIANCE 

An elegant method of introducing longitudinal boost 
invariance with the longitudinal velocity field v z = z/t 
makes use of the notation of general covariant derivatives. 

In an arbitrary coordinate system the equations of mo- 
tion can be written as 



rpTi 



= 0. 



J 



= 0. 



(Al) 



where the semicolon indicates a covariant derivative. For 
tensors of rank 1 and 2 it reads explicitly 



i L = i l + V i i k 
J\p J ,p * x pk J 



Tik rplk 
\p — J- ,p 



pm 



pfc rpim 

pm 



(A2) 
(A3) 



where the komma denotes a simple partial derivative and 
the Christoffel symbols r|- are given by derivatives of the 
metric tensor g ab (x): 



r ij - l9 ks (gikj + 9jk,i - 9ij,h) ■ 



(A4) 



We use this with the following transformation from 
Cartesian to light cone coordinates: 



x^ = (t,x,y,z) 
t = t cosh rj 

z = t smb. rj 



x m = (r,x,y,ri) 

T = 



1 1 t+Z 



(A5) 
(A6) 



In the new coordinate system the velocity field (after 
inserting v z — z/t) is given by 



(A7) 



with Vi = Vi cosh 7/, i = x,y, and 7=1/ ^JY—v^—v 2 . 
Now we turn to the metric of the new system. We have 



ds 2 = g^dx^ dx v — dt 2 — dx 2 — dy 2 — dz 2 

= dr 2 - dx 2 - dy 2 - r 2 dr] 2 

and therefore 
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-T 2 



(A8) 



(A9) 



The only non- vanishing Christoffel symbols are 



r/T 



pJ7 



(A10) 



Finally, by making use of the relations T T% = v{T TT + 
ViP and T m = p/r 2 the energy-momentum conservation 
equations ( |Al| ) turn for n = T,x,y,r] into 



T TT , T + (v x T TT ), x + (v y T Tr ), v 
P+T T7 



(v x T TX ) :X + (v y T T *),y = 

-P.x , 



+ (V X T TV ) tX + {VyT T y) t y 

rpry 

= ~P>y - — > 



while the current conservation (Al) becomes 

f 

f.r + (Vxf),x + (Vyf), V =-—■ 



(Alia) 
(Allb) 
(Allc) 

(Alld) 
(A12) 



We note the explicit appearance of r on the r.h.s. of 
the differential equations, reflecting the dilution of the 
matter due to the boost-variant longitudinal expansion. 
Connected with this is the initial equilib ration time To as 
one of the model parameters. Equation (Alld) expresses 
the fact that, due to longitudinal boost-invariance, the 
evolution is 77-independent. 

Multiplying these equations by r and introducing the 
scaled quantities J p = T7 7 ', T^ v = tT^ v , and p = rp 
leads to the simple form (2.4). 
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